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ABSTRACT 

We report laboratory measurements of the absorption coefficient of solid pa?’a-H 2 , 
within the wavelength range from 1 to 16.7 /im, at high spectral resolution. In addition 
to the narrow rovibrational lines of H 2 which are familiar from gas phase spectroscopy, 
the data manifest double transitions and broad phonon branches that are characteristic 
specifically of hydrogen in the solid phase. These transitions are of interest because 
they provide a spectral signature which is independent of the impurity content of 
the matrix. We have used our data, in combination with a model of the ultraviolet 
absorptions of the H 2 molecule, to construct the dielectric function of solid para- 
H 2 over a broad range of frequencies. Our results will be useful in determining the 
electromagnetic response of small particles of solid hydrogen. The dielectric function 
makes it clear that pure H 2 dust would contribute to IR extinction predominantly 
by scattering starlight, rather than absorbing it, and the characteristic IR absorption 
spectrum of the hydrogen matrix itself will be difficult to observe. 
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1 INTRODUCTION 

As a candidate material for interstellar dust, solid hydro¬ 
gen is not new. But it is not well explored. Shortly after 
being proposed as a possible dust constituent (e.g. Wick- 
ramasinghe and Reddish 1968), it was recognised that the 
pure solid would sublimate rapidly under typical interstellar 
conditions (Field 1969; Greenberg and de Jong 1969), and 
interest in the material soon waned. Subsequent studies of 
“condensed” H 2 have been predominantly concerned with 
the distinct case of H 2 molecules adsorbed onto other mate¬ 
rials — either individual atoms/molecules/ions (e.g. Duley 
1996; Bernstein, Clark and Lynch 2013), or else the surfaces 
or porous matrices of other solids (e.g. Sandford and Alla- 
mandola 1993). Such composites differ from bulk solid H 2 
in many ways, including the spectroscopic properties of the 
hydrogen - e.g the Qi line, which was the main focus of 
Sandford and Allamandola (1993), is entirely absent in the 
bulk (see §2 of this paper) - and consequently we do not 
consider them here. 

Pfenniger, Combes and Martinet (1994) suggested that 
spiral galaxies might contain gas in cold, dense, molecular 
clouds which have escaped direct detection i.e. a form of 
baryonic dark matter. In such clouds the H 2 could reach its 
saturated vapour pressure, allowing pure solid hydrogen to 
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precipitate (Pfenniger and Combes 1994; Wardle and Walker 
1999; Pfenniger 2004). Some aspects of those models are con¬ 
troversial, but not the following point: any gas which is cold 
and dense enough to allow H 2 to precipitate would likely 
have escaped direct detection to date. The study of solid 
H 2 formation sites is, therefore, primarily a theoretical en¬ 
deavour at present. But from those formation sites hydrogen 
grains can be injected into the broader ISM (e.g. by contin¬ 
uous stripping of surface material, or by disruptive events), 
where they may be more readily observed. Furthermore, we 
now know that the thermodynamics of the pure material 
are misleading: in interstellar space, solid H 2 would acquire 
a dense layer of surface charges, and these make hydrogen 
grains much more durable (Walker 2013). Given that it may 
be possible for hydrogen dust to survive, sheer weight of 
numbers, relative to any of the metals, motivates an explo¬ 
ration of what such particles might look like. 

The molecular ion H)j" is one aspect of that appear¬ 
ance. It does not feature in the ionisation chemistry of H 2 
in the gas-phase, but is known to form in the irradiated 
solid (Symons & Woolley 2000; Suter et al 2001; Kumada, 
Tachikawa and Takayanagi 2005; Kumada, Takayanagi and 
Kumagai 2006; Kumagai et al 2007); as such it provides a 
signature of the condensed phase. The vibrational lines of 
this species, and its isotopic variant (HD)^, were charac¬ 
terised by Lin, Gilbert and Walker (2011), who reported a 
striking match to the pattern of strong infrared bands of the 
interstellar medium (ISM). 
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A thorough investigation of the properties of H 2 dust 
is called for, and for that we need to establish the dielectric 
function of pure, solid H 2 . This paper presents a preliminary 
description of that dielectric function. Our model is based on 
our own infrared absorption measurements of the solid in the 
wavenumbeiQ range 600-10,000 cm 1 2 (16.7-1 pm), supple¬ 
mented by published data on the lowest rotational transition 
of the solid, which lies outside the range of our spectrome¬ 
ter. In the case of isolated H 2 molecules the infrared lines 
are electric quadrupole transitions, and consequently the ab¬ 
sorptions we have measured are weak compared to those in 
the UV, which correspond to electric dipole transitions of 
the individual molecules. Although we have obtained no UV 
spectra, the strength of those transitions makes it essential 
to include them in any model of the dielectric properties. 
We have therefore represented the discrete UV transitions 
with Lorentz oscillators and the bound-free continuum with 
a published analytic fit to the measured absorption cross- 
section of the H 2 molecule. We then synthesised the dielec¬ 
tric function via the Kramers-Kronig relations. 

Although our model of the UV absorptions is incom¬ 
plete, the omissions are not very significant in respect of 
the overall dielectric response — the model captures 95% 
of the oscillator strength, for example. However, it is not a 
good description of the detailed structure of the UV absorp¬ 
tion bands of the solid, because it is based on the measured 
response of gas-phase molecules and then simply scaled to 
the density of the solid. This approach fails to account for 
the line broadening which arises from intermolecular interac¬ 
tions within the solid. Away from the UV absorption bands 
themselves, our model should be reliable. In particular it 
should give a satisfactory account of the influence which the 
strong UV bands have on the infrared dielectric function, 
and it thus provides appropriate context for our infrared 
measurements of the solid itself. Despite the shortcomings 
of our dielectric function in the UV, we present the com¬ 
plete model here. The motivation for doing so is simply that 
a better model is not available at present. If it is used in the 
UV, it should be used with caution. 

The structure of this paper is as follows. In the next sec¬ 
tion we describe our laboratory techniques and the infrared 
(IR) absorptions of solid para- H 2 which we have measured]^] 
In section 3 we present our model of the UV response of the 
solid, based on the discrete and bound-free transitions of H 2 
reported in the literature, and this is combined with our IR 
absorption data to yield a single description of the dielectric 
function. Discussion and conclusions follow in §§4,5. 


1 Throughout this paper we employ wavenumbers measured in 
cm -1 , i.e. just the reciprocal of the wavelength measured in cm. 

2 The ortho/para designation refers to odd/even angular momen¬ 
tum quantum numbers. Because ortho-para state transitions are 
forbidden, conversion between the two forms is slow, and these 
sequences constitute separate populations in statistical equilib¬ 
rium. However, at low temperatures all the molecules eventually 
end up in the J = 0 state, so it is the para form which is relevant 
here. Henceforth we use the shorthand p- H 2 for para- H 2 . 


2 LABORATORY MEASUREMENTS OF 

INFRARED ABSORPTION 

At low temperatures, infrared absorption spectra of p- H 2 in 
the gas phase are simple, consisting of the pure rotational 
line So(0), at approximately 354cm _1 (28.2 pm), and a se¬ 
quence of rovibrational lines S„(0), with v = 1,2,3 .... Here, 
and subsequently, we follow the usual spectroscopic nomen¬ 
clature in which changes in the angular momentum quantum 
number by amounts A J = —1,0,1,2,4 are denoted by the 
letters P, Q, R, S, U, respectively. The subscript on the des¬ 
ignation indicates the vibrational quantum number of the 
upper state, and the number in parentheses denotes the ro¬ 
tational quantum number of the lower state, which is always 
0 for P-H 2 at low temperatures. Absorptions to vibrational 
states higher than v = 1 arise from anharmonicity in the 
interatomic potential and are comparatively weak, with os¬ 
cillator strengths that decrease by a factor ~ 10 for each of 
the first few steps up the vibrational ladder (Turner, Kirby- 
Docken and Dalgarno 1977). 

Once P-H 2 molecules are condensed into the solid phase, 
however, much more structure is seen (e.g. Gush et al 
1960; Buontempo et al 1982; Mengel, Winnewisser and Win- 
newisser 1998; Mishra et al 2004). For example, a Ui(0) line 
appears with significant intensity, induced by intermolec¬ 
ular interactions of a central P-H 2 molecule with the sur¬ 
rounding P-H 2 matrijj^] (Steinhoff et al 1994). In addition 
to transitions associated with a single molecule, one sees 
double transitions, e.g. Qi(0) 4- So(0), in which two neigh¬ 
bouring molecules make a simultaneous change in quan¬ 
tum state. Triple transitions have also been identified (Men¬ 
gel, Winnewisser and Winnewisser 1998) but they are much 
less prominent. And each transition manifests a broad side¬ 
band, on the short wavelength side, in which phonons carry 
off energy into the lattice. Phonon-assisted transitions are 
designated with the subscript R. Some transitions which 
are strictly forbidden in isolated P-H 2 molecules can, with 
phonon assistance, acquire substantial intensity in the solid. 
A good example is the Qr( 0) band, which is prominent in 
solid P-H 2 , despite the fact that J = 0 —>■ J = 0 is a for¬ 
bidden transition (Field, Somerville and Dressier 1966); and 
the Qi(0) line itself is absent even in the solid (Hinde 2003; 
Anderson et al 2002). 

A theoretical basis for describing the infrared specrum 
of solid P-H 2 was pioneered by Van Kranendonk and col¬ 
laborators (e.g. Van Kranendonk and Karl 1968). Many of 
the observed features can be understood within a framework 
where the quanta of molecular rotation and vibration are to 
some degree delocalised in the lattice, and the terms “roton” 
and “vibron” are applied to those lattice excitations. 

Those lines which correspond to transitions of individ¬ 
ual molecules are seen to occur at slightly different wave¬ 
lengths in the solid — a “matrix shift”. For example in solid 
P-H 2 the Si(0) line appears at approximately 4,486 cm -1 
(2.229 pm), roughly 12 cnY 1 lower in wavenumber than the 
same line in the gas phase. The IR transitions of solid p- 
Ho are typically much more intense than the analogous gas 
phase transitions, which occur by quadrupole radiation, be¬ 
cause the transitions acquire some electric dipole character 

3 The neighbouring molecules each have a non-zero electric 
quadrupole moment. 
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via intermolecular interactions in pairs of H 2 molecules — 
analogous to collision-induced IR activity (Raston and An¬ 
derson 2007). 

2.1 Experimental method 

We prepared P-H 2 crystals using the rapid vapor deposition 
(RVD) method of Fajardo and Tam (1998). In this method 
the crystal is grown by deposition of pre-cooled P-H 2 gas 
onto a 2 mm thick BaF 2 optical substrate held at ~ 2.5 K 
within a sample-in-vacuum liquid-He cryostat. In order to 
minimise the effects of high-order internal reflections in the 
optical substrate, its front and back surfaces are not paral¬ 
lel; the wedge angle is 1°. For the same reason, we probe 
the sample with a weakly converging IR beam rather than a 
collimated beam. The P-H 2 gas is enriched from normal- H 2 
gas to greater than 99.97% P-H 2 using a variable tempera¬ 
ture ortho/para converter operated near 14.0 K. The ortho- 
H 2 concentration in the sample can be determined using 
the Qi(l) or Ui(l) transitions and their known absorption 
strengths. 

The vacuum chamber that surrounds the cold optical 
substrate has a nominal pressure of 10~ 6 Torr prior to any 
cryogens being introduced. This is sufficiently large that all 
samples become increasingly contaminated with H 2 O and 
CO 2 impurities over the lifetime of the crystal (approxi¬ 
mately 4 to 10 hours). These trace impurities do not cause 
any great difficulty for us, however, as most of their absorp¬ 
tion bands do not overlap with those of the P-H 2 itself. 

High-resolution Fourier Transform IR (FTIR) spec¬ 
troscopy was performed on each sample using a normal- 
incidence transmission optical setup. The IR beam is fo¬ 
cused through the sample with off-axis parabolic mirrors of 
8” effective focal length. The optical path outside the cryo¬ 
stat and spectrometer is contained within an optical conduit 
purged with dry air to minimize atmospheric absorptions. 
FTIR spectra were recorded at various resolutions ranging 
from 0.01 to 0.1 cm -1 both during and after deposition. The 
FTIR spectrometer (Bruker IFS 120HR) used to record IR 
spectra is equipped with a variety of sources, beamsplitters 
and detectors to cover the full range of the spectrum. Mid-IR 
spectra are recorded with a glowbar source, a Ge-coated KBr 
beamsplitter, and a liquid nitrogen cooled HgCdTe detector 
in the region from 600 to 6,000 cm -1 at 0.1 cm -1 resolu¬ 
tion, and with an InSb detector from 2,000 to 5,000 cm -1 
at 0.09 cm -1 resolution. Near-IR spectra are recorded with a 
tungsten source, a CaF 2 beamsplitter, and a liquid nitrogen 
cooled InSb detector in the region from 2,000 to 10,000 cm -1 
at 0.09 cnU 1 resolution. 

The three instrumental configurations just described 
provide high resolution spectra over a broad wavelength 
range, including the fundamental and first overtone vibra¬ 
tional modes. Unfortunately, our spectral coverage does not 
include the lowest pure rotational transition, i.e. So(0), nor 
the peak of the associated phonon branch Sr(0) 

2.2 Thickness determination 

In their original paper in 2001, Tam and Fajardo measured 
the integrated absorption intensities of the Qr(0) + So(0) 
(~4,508 cm -1 ) and Si(0) + So(0) (~4,840 cm -1 ) double 



Figure 1. Upper panel: integrated intensity of the Qi(0) + So(0) 
transition, versus that of the Si(0) + So(0) transition of solid p- 
H 2 , as measured with the three beam-splitter/detector combi¬ 
nations used in the present study. Blue circles denote measure¬ 
ments made with CaF 2 /InSb; green triangles for KBr/InSb; and 
red squares for KBr/HgCdTe. The plotted line is the best-fit re¬ 
lationship, in a least-squares sense, of a proportionality between 
the two integrated intensities. Lower panel: residuals relative to 
the best-ft relationship. 

transitions of p- H 2 and calibrated them against the sample 
thickness determined interferometrically. They also showed 
that the integrated intensities of these double transitions are 
insensitive to crystal structure and the presence of impu¬ 
rities. Furthermore the Qi(0) + So(0) transition is approx¬ 
imately ten times stronger than the Si (0) + So(0) transi¬ 
tion, so one can use the stronger transition for thin sam¬ 
ples (d < 0.02 cm), the weaker one for thick samples (up to 
10 cm), and the combination for intermediate thicknesses, 
allowing for consistency checks. We have therefore followed 
Tam and Fajardo (2001) and used these double transitions 
to determine the thickness of our samples. 

To check the internal consistency of this method, we 
examined the correlation between the integrated intensities 
of these double transitions determined using all three of our 
beam-splitter/detector combinations; results are shown in 
Figure 1. Both FTIR setups that use the KBr beam-splitter 
showed a slope consistent with the previous measurements 
of Tam and Fajardo (2001). However, for the setup that em¬ 
ployed a CaFi beam-splitter, the determined slope was ap¬ 
proximately 6% larger than the literature value and outside 
the 95% confidence range. The reasons for this discrepancy 
are unclear. However, to determine the thickness using a 
common protocol over the entire spectral range we decided 
to use the integrated intensity of the Qi(0) + So(0) double 
transition of 90(±2) cm -2 (Fajardo 2011), and to scale the 
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Figure 2. Line profiles for the Si(0) + So(0) double transition 
(left panel), and the Ui(0) single transition (right panel), of solid 
p- H 2 . These profiles are taken from spectra of the same sam¬ 
ple (approx 3 mm thick), recorded at different times during the 
experiment. Trace (a) was recorded at 1.7 K, in the freshly de¬ 
posited sample; trace (b) during annealing of the sample, for 30 
minutes at 4.3 K; and trace (c) after annealing and re-cooling the 
sample to 1.6 K. The numerical value adjacent to each trace is 
the integrated intensity of each line (in cm' 1 ). The double tran¬ 
sition exhibits some profile evolution, but its integrated intensity 
is steady to within ~ 1%. By contrast the integrated intensity of 
the single transition increases substantially during annealing, as 
fee crystal domains are converted to the more thermodynamically 
stable hep structure. 


integrated intensity of the Si(0) + So(0) accordingly for the 
three setups. 

2.3 Annealing of samples 

Samples of P-H 2 prepared by RVD contain a mixture of face 
centered cubic (fee) and hexagonal close packed (hep) struc¬ 
tures in the freshly deposited material. The fee crystal do¬ 
mains are, however, only metastable: they arise because the 
rapid deposition of material does not allow time for relax¬ 
ation into the thermodynamically more stable hep structure. 
Although a freshly deposited solid p- H 2 sample can contain 
upwards of 50% in fee domains, annealing of the sample at 

4.3 K for a period of several minutes will almost completely 
convert the fee into hep structures. In an astrophysical con¬ 
text there is plenty of time for solid P-H 2 to relax into its 
minimum energy configuration, so it is the more stable, hep 
structure that is relevant, and consequently all of our sam¬ 
ples underwent annealing. Spectra were recorded during de¬ 
position and annealing, and after the annealing cycle was 
complete. 

Although the fee and hep structures of solid p- H 2 are 
very similar, their spectra are not identical. A key difference 
is that in the fee structure, but not in the hep structure, each 
P-H 2 molecule occupies a crystal site with center-of-inversion 
symmetry, so that the induced dipole moment for the p- H 2 
molecule in this site cancels out to a large degree for the 
fee structure. Consequently some of the single transitions 
of solid P-H 2 are much stronger in hep than in fee; a good 
example of this is the Ui(0) transition, which is only induced 
in hep crystal sites (Raston et al 2013). 

Considerations of the symmetry of an individual site in 


the lattice are less relevant to the strength of double tran¬ 
sitions, and those transitions may exhibit integrated inten¬ 
sities which are independent of the fcc/hcp content of the 
sample. That is one reason why double transitions are pre¬ 
ferred for thickness determination — as described in §2.2. 
Figure 2 illustrates both the Si(0) + So(0) and the Ui(0) 
transitions, recorded simultaneously in a single sample: right 
after deposition was complete (trace a), while annealing the 
sample at 4.3 K for 30 minutes (trace b), and after re-cooling 
the sample to 1.57 K (trace c). In detail the profile of the 
Si(0) + So(0) double transition is seen to change slightly 
between the as-deposited, and annealing or annealed spec¬ 
tra; but the total intensity, integrated over the profile, is 
constant to within ±1%. By contrast the Ui(0) transition 
increases in strength during annealing, and is 70% stronger 
in the annealed sample compared to the freshly deposited 
material. 


2.4 Measured absorption coefficients 


From the measured intensity, I{v), we can determine the 
absorption coefficient, a(u), which is a material property 
and thus independent of the sample thickness, d. The two 
are related by 


a{u) 





(1) 


where Io(v) is the measured intensity in the absence of a 
P-H 2 sample. Equation 1 would be exact in the case of a 
material with unit refractive index. As we will see in due 
course, solid p-ffa has a refractive index that is very close 
to one throughout the range of our measurements, so the 
approximation in equation 1 is a good one. For the purpose 
of deriving absorption coefficients, we have used only spectra 
from fully annealed crystals of p- H 2 . 

Figure 3 shows our measured absorption coefficients in 
three distinct spectral regions. The regions around the fun¬ 
damental and first vibrational overtone exhibit the huge 
range of linewidths that is one of the remarkable features 
of the IR spectra of solid hydrogen. For example, the sharp 
transition Ui(0) at 5,261 cm -1 (1.901 pm) has a FWHM 
~ 0.16 cm -1 , whereas the phonon-assisted transition Qb( 0), 
which has a similar peak absorption coefficient (located near 
4,230 cm -1 ), displays a FWHM in excess of 100 cm -1 . 

The Ui(0) excitation is essentially localised on a single 
molecule, because the vibrational energy of this transition 
lies below the P-H 2 vibron band. In turn that is because of 
centrifugal distortion of the H 2 molecule in the J — 4 state. 

The Ui(0) transition has its own phonon side-band 
which is broad and peaks roughly 45 cm -1 higher in 
wavenumber than the zero-phonon line. This phonon band 
overlaps with the double transition Qi(0) + Uo(0), which 
is the source of the sharp peak at 5,320cm -1 (1.880pm). 
Broad, asymmetric, phonon side-bands can be seen to be 
situated in a similar relationship to essentially all of the ro- 
vibrational transitions - both single and double - in figure 
3. 

The double transitions have complicated line profiles 
and may themselves be quite broad (i.e. even the zero- 
phonon transition), as can be seen with the Si(0) + So(0) 
transition in figures 2 and 3. 
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Figure 3. The measured absorption coefficient of solid p- H 2 , 
showing: pure rotational transitions (top panel); the fundamen¬ 
tal vibrational band (middle panel, with inset showing the full 
range of the Si(0) and Qi(0) + So(0) absorptions); and the first 
vibrational overtone (bottom panel, with inset showing the full 
range of the Qi(0) + Si(0) absorption). Assignments are shown 
for the most prominent transitions. All zero-phonon absorption 
fines are accompanied by a broad phonon branch at slightly higher 
wavenumber. Note also the strong Qr( 0) phonon branch, peak¬ 
ing at approximately 4,230 cm ', but no evidence of the forbid¬ 
den Qi(0) line itself (near 4,150 cm -1 ). The continuum absorp¬ 
tion seen below 700 cm -1 is the phonon branch of the So(0) line, 
which lies beyond the A -1 ~ 600cm -1 instrumental limit. No 
wavelet filtering has been applied to the data at this stage. 


It is no surprise that the transitions in the vibrational 
overtone region are generally weaker than those of the fun¬ 
damental band, but some significant peaks in the absorption 
coefficient are nevertheless manifest in the bottom panel of 
figure 3. Anharmonicity of the H 2 molecule is evident in the 
displacement of the Q2(0) + So(0) transition, relative to the 
Qi(0) + Si(0) transition, of more than 200 cm -1 . 

2.5 Data processing 

Although our laboratory spectra provide a wealth of infor¬ 
mation on the absorption coefficient of solid p- H 2 , the data 
in their raw form are not ideal for constructing a model di¬ 
electric function. In this section we describe how we used 
the data to arrive at a model for the infrared properties of 
the solid. 

2.5.1 Saturated data 

Our principal dataset, obtained for the express purpose of 
constructing the dielectric function of solid P-H 2 , consists 
of two spectral scans, covering the range 2,000-10,000 cm -1 , 
of a solid P-H 2 crystal of thickness 3.1 mm. In these data 
we observed that two P-H 2 features - the Si(0) line and the 
sharper peak of the Qi(0) + Si(0) transition - are so strong 
and narrow that they are saturated. We anticipated this dif¬ 
ficulty and therefore recorded similar scans of a crystal of 
thickness 0.31 mm, in which these lines are not saturated. 
Using the spectra of the thinner crystal it was straight¬ 
forward to determine the line profiles in the narrow re¬ 
gions where saturation is a problem in our main dataset. 
Both regions are well represented by Lorentzian profiles 
with FWHM of approximately 0.08 cm -1 . The correspond¬ 
ing model profiles were used to replace the saturated points 
in our main dataset. 

2.5.2 Contaminants 

As noted earlier, our samples of solid p- H 2 exhibit a low 
level of contamination, notably from H 2 O and CO 2 , and 
these impurities contribute to the absorption spectra. How¬ 
ever, the spectral features of these molecules are well known 
and thus easily distinguished from those of P-H 2 . Mostly 
their contributions lie in regions where the p- H 2 absorp¬ 
tion is in any case too small to be measured in our ex¬ 
periments, and in those regions we simply do not use the 
data so the contamination is irrelevant. The exception is a 
combination line of water at 5,338.54 cm - , which is super¬ 
imposed on the phonon continuum associated with the Ui(0) 
and Qi(0) + Uo(0) transitions. This water line is well repro¬ 
duced by a Lorentz oscillator, and we therefore removed the 
feature by subtracting the best-fit Lorentzian profile. 

In addition to the low-level contamination of our crys¬ 
tals, just described, the spectrometer itself - specifically, the 
HgCdTe detector - has, over the years, acquired a surface 
contaminant. The nature of this contamination is unknown, 
but the effect on our acquired spectra is to create strong 
non-linearities, and thus large artefacts in the low-frequency 
absorption coefficients. In practice this meant that the low- 
frequency spectra we acquired specifically to measure the ab¬ 
sorption coefficient were not useful. Instead we made use of 
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some archived spectra from a study of acetonitrile (CH 3 CN), 
at a concentration of 23 ppm, in a solid p-Ho matrix. These 
spectra were taken before significant detector contamination 
had built-up. Strong, narrow (~ 0.1cm -1 ) acetonitrile ab¬ 
sorptions are present in these spectra — e.g. at 919 cm -1 
and 1041cm -1 . But from the known gas-phase spectra of 
this molecule (Nishio et al 1995), the only contaminating 
line expected within the portion of data which we make use 
of is an overtone line at approximately 718 cm -1 . There is 
no sign of a narrow absorption feature at this location in 
our data, so we conclude that it is at or below the level of 
the noise in the data, and therefore insignificant. 

2.5.3 Wavelet filtering 

Our absorption coefficient spectra, as recorded, include some 
low-amplitude, high frequency fringes — evident around 
750 cm -1 in figure 3. These fringes are very likely a residual 
systematic calibration error arising from multiple internal 
reflections in the BaF 2 substrate. Our spectra, particularly 
the top and bottom panels of figure 3, also exhibit significant 
measurement noise. Where there is strong absorption the 
measurement noise is unimportant, but in spectral regions 
where the absorption is weak the signal we’re interested in 
may be masked by the noise. Furthermore, the presence of 
noise causes numerical difficulties in constructing the real 
part of the refractive index, via the Kramers-Kronig rela¬ 
tions, because one expends a great deal of time carefully 
integrating all the details of the noise in the vicinity of the 
point of interest. 

To address these problems we “denoised” our absorp¬ 
tion coefficient spectra using wavelet transforms. Wavelet 
filtering is preferred over spectral averaging because the 
wavelet approach preserves the profile of strong, narrow fea¬ 
tures in the data, whereas spectral averaging does not. 

Within the Mathematica software package, we utilised 
the discrete wavelet transform, with a Haar wavelet; we ze¬ 
roed all coefficients below a suitably chosen threshold, and 
then applied the inverse discrete wavelet transform. The 
thresholding value was chosen to be close to the level of the 
noise in the data, so as to filter out the noise while preserv¬ 
ing the signal in so far as possible. In regions where the ab¬ 
sorption is small, the filtered spectra are piecewise-constant 
approximations to the true, smoothly varying coefficient. 

2.5.4 Model phonon band 

In regions where the absorption coefficient of solid P-H 2 is 
very small indeed, even our “denoised” data do not yield 
an accurate value. In this case there is little merit in using 
the data themselves, and instead we have employed a simple 
model of the continuum absorption based on the observed 
properties of the phonon side-bands. 

In our data Qr( 0) is the best characterised of the var¬ 
ious phonon side-bands: it has a high peak value, and is 
the dominant source of absorption over a range of approxi¬ 
mately 300 cm -1 around the peak. We have therefore used 
this band as a template for describing other phonon bands 
in the dataset. Specifically, the Qr( 0) band provides us with 
a description of the shape of a phonon band absorption; to 
match a particular phonon band we then scale the template 


according to the value of the peak absorption coefficient, 
and the location of that peak relative to the origin of the 
band (i.e. the zero-phonon line). 

The phonon side-band shapes actually differ substan¬ 
tially amongst the various transitions in our data, and our 
procedure is just a simple heuristic which allows us to ex¬ 
trapolate from regions where the absorption coefficient is 
well-measured into regions where it is not. 

In fact the Qr( 0) band itself required extrapolation in 
order to provide a template that is valid over a suitably 
large range of wavenumber. Our Qr( 0) data display approx¬ 
imately exponential decline in the absorption coefficient at 
large distances from the peak of the band, and we there¬ 
fore used exponential functions to extend our phonon band 
template to wavenumbers smaller and larger than those for 
which we were able to measure the band shape. 

2.5.5 Data segments utilised 

The infrared absorptions of solid p-H -2 are clustered around 
the various S„(0) lines (w = 0 , 1 , 2 ...), as is clear from figure 
1 of Mengel, Winnewisser and Winnewisser (1998). Conse¬ 
quently we only utilised our data to define the imaginary 
part of the refractive index in those segments where the 
absorption signal is strong. Those segments are the follow¬ 
ing: 689-850 cm -1 , 4,000-5,450 cm -1 , 8,395-8,790 cm -1 , and 
8,970-8,974 cm -1 . The data for each segment can be seen 
in figure 3. Outside of those segments (but within our in¬ 
strumental window of 600-10,000 cm -1 ), our data were used 
indirectly. That is, the data were used to dehne models of the 
absorption spectrum, and we then employed the models in 
our refractive index description. The models are Lorentzian 
profiles, in the case of the Uo(0) and S 2 ( 0 ) lines - both of 
which are very narrow - and phonon band absorptions as 
described in the previous section. 


2.6 Long wavelength absorptions 

The pure rotational line, So(0), and the peak of the as¬ 
sociated phonon branch, Sr( 0), lie below the 600 cm -1 
wavenumber limit of our spectrometer and are thus absent 
from our data. These transitions are known to be strong ab¬ 
sorptions of solid P-H 2 , relative to the IR transitions we have 
measured, and it is therefore desirable to include a represen¬ 
tation of these features in our description of the dielectric 
properties. 

Our data show the Si(0) line to be well represented by 
a Lorentzian absorption profile with FWHM of 0.084 cm -1 , 
and we therefore represent the So(0) transition with a 
Lorentz Oscillator model (as detailed in §3.1), yielding the 
same profile as the Si(0) line. The wavenumber at resonance 
we use for the So(0) line is 355.7cm -1 (Buontempo et al 
1982), and the integrated line-strength (Buontempo et al 
1982; Okumura, Chan and Oka 1989), taken together with 
our profile shape, indicates a peak absorption coefficient of 
159 cm -1 . 

Unlike the zero-phonon So(0) line, the phonon branch, 
Sr(0), is not expected to follow a simple analytic form. We 
have therefore chosen to represent the Sr(0) absorption pro¬ 
file by scaling the Qr(0) phonon band, as described in §2.5.4. 
Laboratory measurements of the Sr(0) branch indicate that 
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the peak absorption is 5.3 cm -1 , and that the peak is approx¬ 
imately 55 cm -1 above the So(0) line centroid (Buontempo 
et al 1982), and we scaled our template accordingly. 


2.7 Short wavelength absorptions 

There are, of course, rovibrational absorption bands of solid- 
p- H 2 which lie at wavenumbers higher than we can mea¬ 
sure with our spectrometer, corresponding to the second, 
third, fourth etc. vibrational overtones. As already noted, 
those transitions are relatively weak (Mengel, Winnewisser 
and Winnewisser 2013), and we make no attempt to include 
them in our dielectric model. 


3 MODEL DIELECTRIC FUNCTION 

The strongest electromagnetic response of H 2 is to be found 
in the far-ultraviolet (FUV), where various absorption bands 
and the ionisation edge are all located. These FUV transi¬ 
tions are the dominant contributions to the dielectric func¬ 
tion of solid P-H 2 . Measurements of the FUV properties of 
solid H 2 were reported by Inoue, Kanzaki and Suga (1979), 
but there are several reasons why those results are difficult 
to use for our purpose here: they are only presented as op¬ 
tical densities in graphical form; the sample thickness was 
not determined; and the ortho/para content of the samples 
is unspecified. Consequently we are obliged to construct a 
model of the FUV absorptions based on individual molecules 
of H 2 , and scaled to the density of the solid. We first treat 
the contributions of the discrete transitions, and then deal 
separately with the bound-free absorptions. 


3.1 FUV absorption lines 


The lowest-lying absorption bands are the Lyman and 
Werner systems; these are transitions from the electronic 
ground-state to the excited electronic states B(f E+) and 
C( 1 n„) respectively (Field, Somerville and Dressier 1966). 
At the low temperatures of interest to us, essentially all of 
the H 2 molecules are expected to be in the ground rovi¬ 
brational level of the ground electronic state, and the only 
relevant transitions are those which connect to this state. 
This limits us to the R(0)-branch transitions of each ab¬ 
sorption banc 0 i.e. the rotational quantum number of the 
upper state must be J = 1 - but there is no restriction on 
the vibrational quantum number for the upper level. Each 
of these lines can be represented by a Lorentz oscillator, so 
their combined contribution to the susceptibility for radia¬ 
tion of angular frequency uj can be modelled as 


XioM 


ne 2 fj 

eo m e ^ — uj 2 — i uo'jj ’ 

3 J 


( 2 ) 


(e.g. Fox 2012) where the subscript j identifies the contribut¬ 
ing transition, for which fj is the oscillator strength, ujj the 
angular frequency at resonance, and 7 , the damping con¬ 
stant. As usual, e, m e and eo denote the charge and mass of 


4 In contrast to the IR transitions of H 2 , which are quadrupolar, 
the dominant UV transitions are electric dipole, so AJ = 1. 


the electron, and the permittivity of free-space. The quan¬ 
tity n is the number of molecules per unit volume. 

To make use of equation (2) we need values of 
{cjj, f -j , 7 , }. We have compiled suitable data from the liter¬ 
ature, and in Appendix A we list values of these quantities 
for each of the transitions included in our model. In sum¬ 
mary we have 99 transitions, made up of: 38 transitions in 
the Lyman Band; 14 transitions in the Werner Band; and 
47 transitions involving a further 5 excited electronic states. 

The combined oscillator strength of all transitions 
should equal the number of electrons in each molecule, i.e. 
2 , whereas the discrete transitions in our model contribute 
y fj — 0.829. Certainly our list of discrete transitions is 
incomplete, but the major part of the discrepancy is due to 
bound-free transitions. Because they form a spectral contin¬ 
uum, these transitions are not naturally described in terms 
of Lorentz oscillators. 


3.2 Bound-Free continuum 


To model the bound-free continuum we use a piecewise ana¬ 
lytic fit (Yan, Sadeghpour and Dalgarno 1998), to the mea¬ 
sured bound-free cross-section of the H 2 molecule, atf(uj) 
(Samson and Haddad 1994). Near the ionisation edge the fit 
we employ is the modified one due to Liu and Shemansky 
(2004). For low number-densities of molecules, the imaginary 
part of the bound-free contribution to the refractive index, 
is then given by Im{n b ^( c<j)} = nabf(uj)c/2uj. The 
corresponding real part can be obtained from the imaginary 
part by applying the appropriate Kramers-Kronig relation 
— as per Draine and Lee (1984), for example. Knowing both 
real and imaginary parts of the refractive index contribution, 
we can determine the bound-free susceptibility contribution 
from x bf = n If ~ 1 - 

The total electromagnetic response of a dilute collection 
of H 2 molecules, due to both line and continuum interac¬ 
tions, is then just 


x d H = XioH + Xb f (v)- 


(3) 


We have not yet specified the number density of molecules, 
n; that is fixed by the mass density of the solid 
(0.087gcnD 4 ). In gas-phase the dielectric constant scales 
linearly with n, so that e = 1 + x d ■ But in the solid phase 
the packing is so close that the polarisation of neighbouring 
molecules has a significant effect on the electric held at any 
given lattice site, and x = e — 1 7 ^ Xd ■ We model this in the 
usual way, utilising the Clausius-Mossotti/Lorenz-Lorentz 
relation, which yields 


e(w) = 1 + 


3XdH 

3 -XdH 


(4) 


for the solid. Strictly the Clausius-Mossotti/Lorenz-Lorentz 
relation is only valid for cubic crystals, but the preferred 
hep crystal of solid P-H 2 is, as noted earlier, very similar 
to a fee structure, so we do not expect to incur a large er¬ 
ror in utilising the relationship here. (A rough estimate of 
the magnitude of the error can be obtained by analogy with 
solid 4 Hc, which also forms an hep crystal at low temper¬ 
atures. For that crystal the optical birefringence (Vos et al 
1967; Kronig and Sonnen 1958), indicates a departure from 
equation (4) that is ~ 10 -4 of e — 1. This is much smaller 
than the systematic errors involved in our estimate of Xd-) 
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Table 1. The model dielectric constant of solid P-H 2 as a function 
of wavenumber, running over wavelengths from one metre down 
to slightly below one Angstrom. Only the first-two and last-two 
lines are shown here; the full table is available in electronic form. 


1/A 

Re{e — 1} 

Im{e — 1} 

(cm 7 ) 



1.00000 x 10 -2 

0.269283 

1.00205 x 10“ 13 

1.10975 x 10 -2 

0.269283 

1.11202 x 10 -13 


9.93305 x 10 7 

-4.49018 x lO" 7 

2.69114 x 

10 -13 

1.01644 x 10 s 

-4.28809 x 10~ 7 

2.42858 x 

10- 13 



3.3 Infrared absorptions 

Our model of the dielectric function expressed by equations 
3 and 4 does not include the contribution due to infrared 
transitions. As we have measured the IR absorptions (§2), 
we can determine the contribution of the IR transitions to 
e in a manner closely analogous to our determination of Xi,j 
in §3.2. In particular we have 

M-U = g (5) 

(where A is the vacuum wavelength), then Re{n ir } is de¬ 
termined from Im{m r } by applying the relevant Kramers- 
Kronig transform, and the susceptibility contribution is 

Xir ^ir 1 - 

Because we determined the IR absorption coefficient of 
the solid itself, rather than a dilute collection of molecules, 
it is not necessary to apply a local-field correction to the IR 
susceptibility. Our model dielectric constant therefore fol¬ 
lows from adding \ iT to equation 4: 

e(w) = 1 + — *-d( ) , + (6) 

3 -X d {u) 

where x,i is again given by equation 3. The resulting re¬ 
fractive index is shown in figure 4. The same information is 
given in table 1 as e — 1 versus A -1 ; this table is very large 
and is provided in electronic form. 

3.4 Comparison of model with data 

Some quantitative checks on the accuracy of our model are 
available. First, the (real part of the) dielectric constant of 
solid P-H 2 at low frequency (kHz regime) has been measured 
to be e(0) — 1 = 0.2833 (Constable, Clark and Gaines 1975), 
whereas our model predicts 0.2693. Secondly, measurements 
of the (real part of the) refractive index of solid P-H 2 in 
the visible and near-infrared bands yield values in the range 
0-13 < \ft— 1 < 0.14 (Perera et al 2011), whereas our model 
yields 0.127-0.132 across the same region of the spectrum. 
Thirdly, in the hard X-ray region the binding energy of the 
electrons is negligible and the refractive index should be just 
that of a cold plasma of electron density 2 n, whereas our 
model tends towards a cold plasma with density 1.894n. This 
is another way of saying that we have captured almost all of 
the oscillator strength. On all three of these tests our model 
performs well; it systematically under-predicts the response 
of the medium at all frequencies, but only by about 5%. 




Figure 4. Real (black) and imaginary (blue) parts of the refrac¬ 
tive index (minus 1) of solid P-H 2 , as a function of wavenum¬ 
ber (top panel). Our laboratory data lie within the region 1.2 < 
— log j q A(//m) 0, and absorptions lying outside this region have 

been modelled with the aid of published data. At low wavenum¬ 
bers the real part of y/e— 1 is positive, whereas at high wavenum¬ 
bers it is negative, and in the absorption bands around 10 ;im _1 
it changes sign many times; we have therefore plotted the absolute 
value. Between 11 and 14.5 eV (0.95 < — log 10 A(/rm) < 1.07) dis¬ 
crete transitions, detailed in Appendix A, are the main contribu¬ 
tion to absorption. The bound-free continuum dominates above 
the ~15.4 eV (—log 10 A(/rm) ~ 1.09) ionisation edge. The red 
dashed line shows a rough confidence level for Im{n — 1}: below 
this line the absorption coefficient is not well characterised. The 
middle and bottom panels show the pure rotational and the fun¬ 
damental vibrational band absorptions, respectively, on a linear 
scale. 
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In comparison with the data of Inoue, Kanzaki and Suga 
(1979), our model for the UV response of solid p- H 2 displays 
discrete transitions which are much too sharp. That is be¬ 
cause our model for the UV is based on a dilute collection 
of H 2 molecules, but scaled to the density of the solid. As 
such the model ignores the intermolecular interactions, oc- 
curing in the solid, that act to broaden the transitions. We 
reiterate, therefore, that our model of the UV resonances 
should not be relied on within the UV band itself - it is 
a good approximation only at much lower and much higher 
frequencies. 


4 DISCUSSION 

Compared to the dielectric constants of conventional in¬ 
terstellar dust materials - i.e. “astronomical silicate” and 
graphite (e.g. Draine and Lee 1984) - the most striking as¬ 
pect of figure 4 is the weakness of the absorption at energies 
below the Lyman band. Because most stars put out most 
of their power in that region, it is immediately clear that 
any pure hydrogen dust in the diffuse interstellar medium 
will be subject to comparatively little radiative heating. At 
present it is unclear how much the absorption evident in 
figure 4 contributes to heating of the crystal, as the elemen¬ 
tary molecular rovibrational excitations are long-lived and 
de-excitation could be principally radiative. Certainly the 
phonon branches contribute heat immediately on absorp¬ 
tion of a photon, but the phonons themselves carry off an 
average < 10 2 cm -1 per absorption, and that is small com¬ 
pared to the energy of the absorbed photon for most of the 
spectral range in figure 4. 

At low wavenumbers the real part of the refractive in¬ 
dex of solid H 2 , although large compared to the imaginary 
part, is also much smaller than that of conventional grain 
materials (e.g. Draine and Lee 1984). This means that the 
wavelength of peak extinction for a pure hydrogen grain is 
expected to be shorter than that for a silicate or graphite 
grain of the same size. 

In addition to the possibility of pure solid H 2 grains, our 
results will be useful for calculating the optical response of 
composite grains consisting, for example, of thick H 2 man¬ 
tles on cores of other materials - refractory or icy. Although 
the H 2 mantles themselves would not absorb much starlight, 
the core material may absorb quite strongly, thus heating 
the whole grain. The radiative heating of such grains in the 
diffuse ISM may thus be too high for an H 2 mantle to per¬ 
sist, and consequently grains of this type are likely to be of 
interest primarily in the context of dense interstellar clouds. 

One of our motivations for precise characterisation of 
the infrared dielectric constant of solid p-IU was the possi¬ 
bility of identifying a spectral signature by which this mate¬ 
rial might be recognised in interstellar dust. In principle we 
have achieved that, as the spectral shape of the absorption 
coefficient in figure 3 is really quite distinctive. However, 
these absorptions are very weak and the imaginary part of 
the refractive index is small compared to the real part, at 
all wavenumbers below the Lyman band. That means that 
it will be difficult to use these IR features to test whether 
interstellar dust contains solid H 2 , because the extinction 
curve should be dominated by scattering and is thus ex¬ 
pected to be relatively featureless (upper panel in figure 4). 


For example: the strongest line in our data is Si(0), and in 
the vicinity of that line the real part of n — 1 varies by less 
than 2% of its local value. Furthermore that fluctuation is 
confined to within ±0.1 cm -1 around the line, so even the 
strongest IR absorptions would be difficult to detect in the 
extinction curve. 


5 CONCLUSIONS 

We have characterised the dielectric constant of solid para- 
H 2 . Our description is based on laboratory measurements 
of the absorption coefficient of the solid in the wavelength 
range 1 < A(/rm) < 16.6, together with published data on 
the low-frequency absorption spectrum of the solid and pub¬ 
lished data on the UV absorptions of the H 2 molecule. The 
real part of the dielectric constant is determined from the 
imaginary part by using the Kramers-Kronig relations. Be¬ 
cause the UV transitions had to be modelled on the basis 
of absorptions measured in gas phase, the spectral structure 
in the vicinity of the UV resonances is not to be trusted. At 
higher photon energies, where the dielectric properties re¬ 
semble a free-electron plasma, the optical constants appear 
to be good to ~ 5%. Comparable accuracy is demonstrated 
for the real part of the dielectric constant at frequencies 
well below the UV resonances. The accuracy of our own 
measurements of the absorption coefficient is also limited, 
by sample thickness determination, to ~ 5% in regions of 
strong absorption, but is much poorer in regions which are 
highly transparent. 

The infrared absorption spectrum of solid p- H 2 is very 
distinctive, with a great variety of line shapes and widths, 
and exhibits transitions which are unique to hydrogen in the 
condensed phase. However, all infrared absorptions are weak 
and consequently the imaginary part of the refractive index 
is small compared to the real part. We therefore expect that 
solid P-H 2 contributes to extinction primarily via a feature¬ 
less scattering continuum, making it difficult to recognise 
hydrogen dust via the IR absorptions of the pure matrix. 
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APPENDIX A: DISCRETE TRANSITION LIST 

In this appendix we provide information on the data used 
in our Lorentz oscillators model of the discrete transitions 
of H 2 , described in §3 of this paper. We supply the data 
themselves as a list of transitions, in the form of an electronic 
supplement, with the following properties. 

There are 99 lines in the list, each one is from the R(0) 
branch, and for each line there are five pieces of data. The 
supplement is thus a table with 5 columns and 99 rows. The 
columns are: (1) the index number we assigned to identify 
the line; (2) the vibrational quantum number, v, of the upper 
level; (3) the wavenumber of the transition, in cm -1 ; (4) the 
Einstein A 21 coefficient of the transition, in MHz; and (5) the 
total decay rate, T, for the upper level, in MHz. The ordering 
of the lines is in groups according to the excited electronic 
state of the transition, and within each group the ordering 
is by increasing values of v. The groups are as follows: 

Lines 1-38: state R( 1 E„), known as the Lyman Band. 
Wavenumbers and A 21 coefficients taken from Abgrall et al 
(1993a), except for v = 23, 33 and v = 35 — 37 which were 
taken from Abgrall and Roueff (1989). All values of T were 
taken from Abgrall, Roueff and Drira (2000). 

Lines 39-52: state C'( 1 n„), known as the Werner Band. 
Wavenumbers and A 21 coefficients taken from Abgrall et al 
(1993b), except for v = 13 which was taken from Abgrall 
and Roueff (1989). All values of T were taken from Abgrall, 
Roueff and Drira (2000). 

Lines 53-59: state B , ( 1 Ej/). Wavenumbers and A 21 co¬ 
efficients taken from Abgrall et al (1994). All values of T 
were taken from Abgrall, Roueff and Drira (2000). 

Lines 60-62: state D( 1 n u ). Wavenumbers and A 21 coef¬ 
ficients taken from Abgrall et al (1994). All values of F were 
taken from Abgrall, Roueff and Drira (2000). 

Lines 63-77: state B" Ej/). All data taken from 
Glass-Maujean et al (2007). Note that this electronic state 
exhibits a double-well structure; the vibrational quantum 
number given in our line list corresponds to the outer well. 

Lines 78-91: state D'( 1 n u ). All data are from Glass- 
Maujean et al (2008a). Note that the total decay rate for 
the levels v = 3 — 9 was smaller than could be measured 
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Table Al. The first-two and last-two lines in the electronic table 
of transitions used in the Lorentz oscillator model of §3. The 
complete table is available in electronic form. 


Line v 1/A A 21 T 

# (cm^ 1 ) (MHz) (MHz) 


1 0 90242.36 3 1,900 

2 1 91558.74 11 1,700 


98 8 133141 3.1 15,000 

99 9 134023 5.1 330,000 


by Glass-Maujean et al (2008a). We have assumed the value 
10 4 MHz for each of these levels, which is a factor of 12 
smaller than the lowest measured value for this state. 

Lines 92-99: state 5pcr( 1 S(/). All data are from Glass- 
Maujean et al (2008b). Note that we adopt their calculated 
value of A 21 for v = 3, as the observed line is blended, and 
for the width of this line we have used the value they list for 
P( 2), v = 3, which corresponds to the same upper level. 

Table A.l reproduces the first-two and last-two lines in 
the list. The various quantities listed in the table are related 
to the quantities employed in the Lorentz oscillator model 
thus: 



Si 


3 . A 2 m e 
4-7T " 1 ctfh ’ 


(Al) 

(A2) 


where a.j ~ 1/137 is the Fine Structure Constant, and h is 
Planck’s Constant. The damping constant is simply 

ii = ^r. (A3) 
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